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We investigate tlie behavior of the dc electronic current J'*'' in an interacting quantum dot driven 
by two ac local potentials oscillating with a frequency f2o and a phase-lag if. We provide analytical 
functions to describe the fingerprints of the Coulomb interaction in an experimental J'*'^ vs ip 
characteristic curve. We show that the Kondo resonance reduces at low temperatures the frequency 
range for the linear behavior of J'*'^ in Qq to take place and determines the evolution of the dc-current 
as the temperature increases. 
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I. INTRODUCTION 

In the beginning of the new century we are witnessing 
an increasing interest towards dc transport induced by 
pure ac fields. Quantum pumps, where transport is gen- 
erated by applying harmonically time-dependent gates 
oscillating with a phase-lag tp at the walls of semiconduct- 
ing quantum dots, areparadigmatic examples realized in 
the laboratory [i|, 0, @]- On the other hand, the possi- 
bility of exploring the Kondo regime in semiconducting 
and carbon nanotube quantum dots provides a unique 
test system to understand the role of electronic corre- 
lations in quantum transport 4, 5]. The combination 
of ac pumping mechanisms with many-body interactions 
constitutes a challenging avenue of research. On the ex- 
perimental side these studies are likely to be feasible in 
the near future, since although these setups employ very 
slowly oscillating fields, great efforts are currently being 
devoted to increase the range of operational frequencies 
Q. The use of superconducting junctions as ac genera- 
tors seems to be a promising methodology in this direc- 
tion [Q]. 

Since the celebrated proposals of Refs. [3,[g|j the "adi- 
abatic approximations" are at the heart of the theoretical 
work on pumping in quantum dots driven at their walls 
B [l3, llll- Within these approximations the induced 
dc-current J'^'^ is proportional to the pumping frequency 
r^Oj describing the regime where flo << t~^ , r being 
the characteristic time for the electrons to travel through 
the dot. Few theoretical studies have addressed the 
problem of many-body interactions in quantum pumps 
[13 . [l3l . [l4| . Most of the work has been centered in 
"adiabatic approximations" [ij, [3j ^-nd the electronic 
interactions are usually included within the slave boson 
mean-field approximation [l3|, which does not properly 
account for inelastic many-body effects. 

In the present work, we also focus on near-equilibrium 
regimes, where JIq is lower than the Kondo tempera- 



ture Tr-, but we explore the effect of the interactions 
beyond the "adiabatic" regime. To this end we com- 
bine two methods: (i) the treatment of the interactions 
by a second-order self-energy as considered in Ref . jl5| , 
which has been a successful tool to study dc transport in 
the Kondo regime [T^ |l3| , and (ii) the Keldysh Green's 
functions formalism with the Fourier representation of 
Ref. [ig|, which has been used to study models of 
non-interacting quantum pumps at arbitrary frequency 
[l6l [l8| . We provide analytical expressions to identify 
the fingerprints of the interactions in the J"**^ vs (/? char- 
acteristic curve, which is the feature usually explored ex- 
perimentally [11]. We anticipate that the development 
of Kondo resonance imposes limits in the range of fre- 
quencies for the J^^ oc JIq behavior to be observed while 
inelastic scattering induced by the interactions tends to 
restore such a behavior. 

The work is organized as follows. In the next section we 
present the model and technical details on the theoretical 
approach. We present results in section IV. Section V is 
devoted to a summary and conclusions. 



II. THEORETICAL FORMULATION 

A. Model 

We describe the quantum pump in terms of a general- 
ized Anderson impurity model (see sketch of Fig. [T|) : 
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FIG. 1: (Color online) Sketch of the setup. The two external 
ac potentials, and the induced potential at the interacting site 
are indicated with solid and dashed arrows, respectively. 



being ei{t) = 5i-i[E + Fcos(floi + ^)] + 5i^i[E + 
V cos{Q,Qt)\. The dot with Coulomb interaction U is in- 
serted between two barriers of height E at which two 
ac fields are applied with amplitude V . The degrees 
of freedom of the reservoirs are denoted with ka, being 
a — L, R, which are at equal chemical potential /x and 
temperature T — 1//3. The reservoirs are attached at 
the positions /„ = ±1 of the central structure through 
hopping terms with amplitude Wc- 



except for the fact that the retarded Green's function cor- 
responds in the present case to the interacting system de- 
scribed by the full Hamiltonian ([1]). Instead, the inelas- 
tic contribution accounts for processes in which electrons 
experiment decoherence originated in the many-body in- 
teractions. It reads: 

r"{k,uj) = 2 ^ Re|r(w) A>(a))S^(fc',cj-mo) 



fe' = -c 



A<(w)I]5'(fc',w- 



-kno) 

k'uj 



{k - k')no) 



ia = ±l 



xGl 



o(fc,w - kflo) I 



(4) 



being A^(w) = «/(w) and A>(w) = ~*(1 — /(^)); which 
depend on the Fermi function /(w) = l/(e^^"~'^' + 1). 
We have also introduced the Fourier representation of 
Ref. [lg| in the lesser and bigger self-energies Sq '^ [k, i^), 
which describe the many-body effects due to the Coulomb 
interaction: 



B. Charge currents 

Following the procedure of Refs. [16], the dc-current 
along the structure can be expressed as a sum of an elastic 
and an inelastic component as follows: 



jdc ^ jel _(_ J^ 



E 

k— — oo 



duj 

2^ 



[r'ik,u;)+I"\k,u;)], 



where, adopting units with e = h = 1, and, to simplify, 
omitting explicit reference to the spin index: 

P\k,uj) = 2\j{uj~kQ.Q)- f{Lo) V{uo)T{uj~knQ) 

X E 'a|Si„,i^(fc,C^-mo)P, (3) 

where V{uj) — SvrlwcPX^fe ^("^ ~ ^feo) is the tunneling 
rate from the structure to the reservoirs and Qi^v {k, lo) is 
the fc-th Fourier coefficient of the retarded Green's func- 
tion: 



+ 00 

k 



E e-^-^«* r°°^e— (*-*')g,.(fc,c.). 



III. TREATMENT OF THE COULOMB 
INTERACTION 

A. Time-dependent Hartree-Fock approximation 

While the above expressions are in principle exact, 
the many-body self-energies must be calculated at some 
level of approximation. The lowest order in the interac- 
tion U corresponds to the time-dependent-Hartree-Fock 
(TDHF) approximation. The many-body problem is re- 
duced to consider Hamiltonian JT]) with [/ = and a 
renormalized level given by 



^TDHF 



(t) = ei{t) + S,,oU 



+00 



TDHF 



{k)e 



ikflnt 



(5) 



/c— — oo 



The Fourier components of the local particle density 
nQ^^^{k) must be evaluated self-consistently. This level 
of perturbation theory leads to a vanishing inelastic con- 
tribution I"''{k,uj) and it is not suitable for the descrip- 
tion of Kondo physics. It is, however, interesting to no- 
tice that this simple approximation already introduces a 
non-trivial ingredient in the pumping problem. Namely, 
the effective emergence of an additional pumping center 
in the interacting site (see Fig. [1]). 



The elastic component takes into account processes 
where electrons propagate coherently along the dot struc- 
ture experimenting virtual absorption or emission of 
quanta of frequency fig at the pumping centers. The en- 
suing expression coincides with the non-interacting one. 



B. Second order self-energy approximation 

In order to go beyond the TDHF description we con- 
sider the second order self-energy (SOSE) given by the 



bubble diagram of Ref. \l3\ and generalize the procedure 
of that work to situations with an harmonic dependence 
on time. Concretely, we consider the following lesser and 
bigger components of the self-energy: 

X [G°^>^<(fc2, t - tTG^'^'^ik -h+k,,t- t'), (6) 

where the propagators are GQ^'^{k,uj) = G^Q'^ikjUj) 
with 



G°;>^<(fc,c^) = J2 ^liA^ + k',Lu- k'Qa) 



a.k' 



xA>'<(cj - k'no)[g?,,i^{k',uj- k'rto)]*, 



(7) 



div 



being the non-equilibrium retarded Green's functions 
Qii,{k,uj), the solution of the Dyson's equation corre- 
sponding to the Hamiltonian ([T]) with U = and 



and it is introduced in the Dyson's equation: for the full 
dressed retarded Green's functions Qi^ii{k,u)). The latter 
are evaluated by using the renormalization procedure of 
Ref. 16] . The algorithm is combined with the use of fast 
Fourier transform between the variables t — t' <-> oj. 



C. Low amplitude and lo^v frequencies expansion 

For small pumping amplitudes and frequencies, it is 
possible to derive an approximate expression for the dc 
current that is accurate up to 0{V^) and up to ©(J^q). 
Such a procedure would correspond to a generalized "adi- 
abatic approximation" within the present formalism and 
will allow us to get analytical expressions to gain insight 
on the behavior of the dc-current. 

As a first step, we truncate the harmonics of 
the induced potentials up to Eg (0) and Eq {1) = 
[^0 (~1)]* — ^eff, and the harmonics of the self-energy 
up to I](0, uj) (higher harmonics involve terms 0{V^) and 
0{U'^V)). The Dyson's equation for the retarded Green's 
function reads: 



er^^{t)=edt)+Si,o E ^o^^(^)^ 



— ikClQt 



(8) 



fc— — C30 



The Fourier components of the effective potential 
Eq (k), are determined self-consistently from the condi- 
tion that the occupation of the interacting site evaluated 
within the SOSE approximation equals the one evaluated 
with (HOl), i.e. 



<^^-(fc) ^ -^ 



duj 

duj 
2^ 



G<{k,uo) 



Gr(fc,c.), (9) 



where the lesser Green's of the first equality contains the 
full dressing by the self-energy: 



Gff{k,uj) ^E^ijjfc-f A',c^-feU) 



(10) 



a..k' 



\><{uj - k'nom'jAk',^;- k'rio)]* + 

Y^ Qifi{k + k' - fc", uj-{k' - k")na) X 

k',k" 

E>'<(fc",c^ ~ k'no)[gi,fl{k',Lj- k'no)]*, 

while that of the second one corresponds to ([7]). The 
above procedure ensures the fulfillment of Friedel- 
Langreth sum rule in the equilibrium limit [15]. Unfortu- 
nately, however, this procedure is not enough to ensure 
the conservation of the inelastic component of the cur- 
rent. 

The retarded self-energy is then obtained from 



G«,(i,c^) ^GlAco) + 



(12) 



J2 Gf:At,u; + no)ef^\l)e-^''«'GlAio) + 

1 
E G5(t,.. - f7o)ef ^■(-l)e^^«*G°,(c.), 

where the pumping potentials contain the external time- 
dependent fields as well as the time-dependent potential 
induced by the interactions: 

ef ^(1) = V[6,,.,e-^^ + (5,, i] + K//<5,.o, (13) 

with [ef ^(1)]* = £^^^(-1). The Green's functions 
G" ;/ (w) correspond to the stationary part of the Hamilto- 
nian. The solution of p^ up to the first order in ej-'{k) 
casts 



Gf,,(i,o.)= Yl Qiri^M^ 



— iujQt 



(14) 



fc=-i 



with 



gi,v{Q,Lu) =Gli,{^) 
1 
g,,,(±l,^) = E G?,,(c.±f]o)£f^(±l)G°,(c^ll5) 



Ys^{k,t-t') ^Q{t-t')[Yj^{k,t-t')-Y.^{k,t~t')], (11) Eq. © leads, within this approximation, to the following 



self-consistency condition to evaluate Veff- 
1 

Veff = U Y. 



duj 



2tt 



ro(a;)/(^)£f^(l) 



+G°,Jc.)G°,,H*G,\(c.-f7o)1. (16) 

A rough estimate of the pumping potential without self- 
consistency is: 



4°)(l)~t/^Ao(l + e-'^), 



(17) 



being 



x[GO,i(c. + f^o)G?^,Jc.)G°,,Jc.)* 
+G^,,J^)G^(c.)*G?,,-J..-f7o)1. (18) 

Thus, within this rough approximation, Veff ~ £o (1) = 
C/yAo[l + cos((/9) — zsin((/3)]. Within the self consistent 
procedure, however, V^f f = V^^r + iV^rr, with both real 

and imaginary parts of the form V^fJ ~ Aq sin((^) + 

A-l^ '{1 + cos{tp)). On the other hand, expanding the 
above expressions up to the first order in Qq we find 

Veff - K°// + KV/"o. 

Expanding ^ in powers of fio i casts for the first order 
contribution: 

jdc _ jei _ f]gr2(/i) 

^ fc{|t;i,i(fc,;,)|2_|(?_i,_i(A:,^)|2 
fc=-i 

+ |ei,-i(fc,M)l'-|e-i.i(fc,M)l'}- (19) 

Substituting the 0th order term in Qq of ([TS]) into P^ . 
it is found 

J^' ^f7or2(;,)^[£f/(l)£f//(l)*- 

ef^(-l)ef^^(-l)1(Ay-A,0, (20) 

being 



A,7 = E ^"G'?^..(A^)GL,i(M)*7/,(^a), 

7^(1) = E G,",Ja*)G°_,Ja^)* (21) 

In the case of the rough estimate for the induced pumping 
potential £q (1) defined above, a relation between the 
current and the phase-lag of the form J*^' ex sin(<p), as 



in the non-interacting case, is obtained. However, when 
the pumping potential is evaluated self-consistently, the 
current behaves as follows: 



r 



r 



no 






V 



Fsin(^)Af 



+ 



(Vo sm{ip){cos{^) + 1) + Vo sin2(^))A^'](22) 

where we have defined F*^' = r(/x), being Fq — 
4|'u;p/F(/i), which is approximately the effective tunnel- 
ing rate from the interacting site to the leads. The dimen- 
sionless functions Xf as well as Vq ' ^ depend on Veff 
and on the parameters Xij of eq. ([2T|) . It is important 
to note that Eq. (|22p reduces to the the non-interacting 
result J'^'^ ^ V'^ sm{ip) originally proposed in Ref. [SJ, 
for V^ = VJj =0. Such a dependence on ip has been ob- 
served experimentally in Ref. [l| and is a manifestation 
of quantum interference between processes of coherent 
emission and absorption of quanta at the two pumping 
centers. In the interacting system the terms oc Xf of 
Eq. (P^ are a consequence of additional interference with 
scattering events at the pumping center induced by the 
interactions. 

The inelastic contribution exhibits the same structure, 
except for a different prefactor F*" = —2Iin[T:^'^(p)], be- 
ing So''(w) - S^(0,a;) the SOSE of the system with 
V" = 0, instead of F'^'. For low T, F™ - (w - i^f, thus, 
for low flo J™ ^ and the elastic processes account for 
the full dc transport. 



IV. RESULTS 

A. Weakly interacting regime 

We begin with the analysis of the effect of the interac- 
tions in the behavior of J^c as a function of JIq at T = 
for a given chemical potential fi and phase-lag tp. In 
Fig. [2] we present the behavior of the different contri- 
butions to J***^ obtained from the numerical solution of 
the full Dyson's equation retaining self-energy compo- 
nents ^^^'^{kjUj) up to |fc| = 2, with the self-consistent 
evaluation of Eg ■' (0), and Sq ■' (1) = Veff- Energies, fre- 
quencies and currents are expressed in units of Fq. The 
inset shows the corresponding local density of states at 
the interacting site pa{uj) — — 2Im[5o,o(0, w)] along with 
the density of states corresponding to the non- interacting 
system {U = 0). The position of the resonant peak ex- 
periments a shift equal to ee//(0). For low U the effects 
of T,{k,uj) are vanishingly small and the description ef- 
fectively reduces to TDHF. The induced Veff is also very 
small {Veff ~ 1 X 10~^ «V). For low C/, the inelastic 
contribution to the current {0{U'^)) is negligible. There- 
fore, within this regime, J'^'^ is qualitatively similar to 
the non-interacting one, shown in dotted line. As the 
hybridization between the dot and the side reservoirs is 
sizable, the resonant peak is wide and the "adiabatic" 
(ex r^o behavior) is observed within a wide range of f2o- 
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FIG. 2: (Color online) r^ (solid line) , J'" (dashed line) 
and J'^'^ (circles) for (^ = 7r/2, as functions of the driving 
frequency fio for [/ = 1. The inelastic component is J'" ~ 0. 
The dc-current for the non-interacting system ([/ = 0) is also 
shown in dotted lines for comparison. The corresponding local 
densities of states at the interacting site po(<^) are indicated 
in the inset (solid lines), along with the non-interacting one 
(dotted lines). Other parameters are _E = 2, w = 2, Wc = 8, 
V = 0.4 and n — 2. (Energies, frequencies and currents are 
expressed in units of Fo ) . 



B. Kondo regime 

Let us now analyze the more subtle Kondo regime, 
which takes place at higher U. In Fig [3] we show a series 
of plots similar to those of Fig. [H but corresponding to 
a value of U for which the Kondo effect takes place. For 
these parameters, the density of states at the interacting 
site po{uj) exhibits the characteristic Kondo resonance at 
uj = with the two high-energy side features centered at 
uj — ±U/2 (see inset). As a function of JIq the behavior 
of J'^'^ significantly departs from the non-interacting one 
in this regime and there are several issues to comment in 
connection to Fig. [31 First, it is clear that self-energy 
effects now play an important role. This is evident in 
the drastic changes experimented by the density of states 
as well as in the fact that the current departs from the 
behavior predicted by the TDHF approximation, which 
is shown in light dot-dashed line in the figure. 

The second feature to remark is that, although inelas- 
tic processes are negligible for low fip, as discussed in 
Section IILC, they become sizable as fio increases even 
at T = (see dashed line of the main frame of Fig. [3]). 

The third remarkable issue is the lose of the "adia- 
batic" behavior, namely, the departure from the linear 
dependence in flo- We identify two ingredients that con- 
tribute to this effect: (i) The first one is the induced 
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FIG. 3: (Color online) J"' (solid line) , /" (dashed line) and 
J"**^ (circles) for ip — 7r/2, as functions of the driving frequency 
flo for U = 10. The details are the same as in Fig. (2] 
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FIG. 4: (Color online) "Transmission function" T{uj) for a 
noninteracting structure with Qq = 0.1 and Wc = 0.64 (upper 
left panel) and Wc — 8 (lower left panel). Other details are 
the same as in Fig. (2] 



pumping potential at the interacting site, which becomes 
sizable Vg// > 1 x 10~^ and changes significantly with 
r^o- This feature manifests itself even at the TDHF level, 
in which case the range of pumping frequencies where 
jdc(TDHF) ^ Q^ becomes very narrow (f^o < 2 x 10"^). 
(ii) The other ingredient is the development of the Kondo 
resonance. Recalling that the underlying assumption for 
low frequency expansions like the one leading to Eq. P^ 
is that the typical width of the energy levels of the struc- 
ture is r^^ >> J7o, it can be understood that the reduc- 
tion of the width of the resonant level due to the Kondo 
effect originates a concomitant reduction of the range of 
r^o for the "adiabatic" behavior to be observed. For the 
parameters of Fig. [31 which correspond to Tr- ~ 0.5, 
such a linear behavior is not captured even close to the 
lowest pumping frequency considered (fio = 1 x 10~^). 
The final issue worth of mention is the inversion in the 
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FIG. 5: (Color online) Integrand /(cj) for Jlo = 0.1 and U = 
10. Different plots correspond to different equidistant values 
of the chemical potential —2 < fi < 0.4. In the inset, the 
mean occupation per spin of the dot is shown as a function 
of the chemical potential. Other parameters are as in Fig. (2] 



sign of the current with respect to the non-interacting 
case. In general, even in the simpler case of a system 
without many-body interactions, the issue of the sign of 
the current is one of the most delicate features to predict 
in a problem of quantum pumping. In the case of struc- 
tures with low hybridization with the contacts presenting 
a landscape of well separated resonances, the issue of the 
sig n of the current has been analyzed in detail in Refs. 
[ifljH)!!^- In those works, the origin for the sign inversion 
of the current has been identified to be the interference 
between two resonant electronic levels mixed by a high 
pumping frequency. In order to gain insight on the dif- 
ferent processes involved in the present problem, let us 
first analyze the behavior of the sign of the current in the 
non-interacting limit of our setup of Fig. [TJ In that case 
the dc current ^ can be also written as follows pjjj : 



j: 



dc 



dujf{uj)T(uj), 



(23) 



where the "transmission" function T{lj) depends on the 
Green's functions of the system with U — Q. Plots of the 
function T{u!) and its integral between — oo and /z, which 
is equivalent to Jnon-int ^^ ^^^'^ temperature for two dif- 
ferent hybridizations Wc, are shown in Fig. 31 In the 
upper panels, corresponding to a small Wc three features 
associated to the three eigenvalues of the tight-binding 
structure with three sites I = —1,0,1, can be distin- 
guished. Instead, for the higher Wc chosen to capture the 
Kondo regime used in Figs. [Hand [31 the three levels of 
the uncoupled structure are mixed and only one feature 
can be distinguished. As a consequence of the ensuing 
combination of quantum states, the transmission func- 
tion as well as the current changes the sign within a wide 
range of // with respect to the weakly coupled structure. 
In the case of the interacting system, we analyze the 
behavior of the function I{uj) — ^^.[/'^'(fc, a;)-t-/'"'(A:,a;)], 
which, when integrated over lo gives the total dc current. 



^ 




FIG. 6: (Color online) J'*'^ for flo — 0.1 as a function of 
the phase- lag ip for [7 = 1 (thick solid line), and U — 10 
(thick dashed line). The current corresponding to the non- 
interacting dot (with [7 = 0) is shown in dotted line, for 
comparison. The symbols are fits with a function Aq am{tp) + 
Aisin^(<^) -I- yl2sin(<^)(l -|- cos{(p)), suggested by eq. (|22 



Inset; Effective potential Veff 
in symbols with the function: 
Other details are as in Fig. [2] 



Vg'jj -I- iV^'ff along with fits 
Aosm{<p) + Ai{l + cos{<p)). 



The behavior of this function within the Kondo regime 
is shown in Fig. [5l Unlike the function T{uj) defined for 
the non-interaction system, the function I{uj) changes as 
the chemical potential changes. For this reason, several 
plots corresponding to different values of the chemical 
potential, for which we have verified that the Kondo res- 
onance is developed, are shown in the figure. For each 
/i there is a feature in /(w) around w ~ /i which pre- 
cisely indicates the electronic transmission through the 
Kondo resonance. Notice that these features are very 
narrow and resemble the lowest energy one of the upper 
left panel of Fig. [H which corresponds to a resonance 
for a dot with low hybridization with the reservoirs. In- 
terestingly, for /i ~ — 1, /(w) experiments a phase shift 
of TT, leading to a concomitant change of sign in the dc 
current. This particular value of the chemical potential 
corresponds to a charge population per spin of the dot 
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'0,cr 



(0) ^ 1/2. 



We show in Fig. [SI the behavior of the dc-current as 
a function of the phase-lag. The limiting case of C/ = 
is plotted in dotted lines in Fig. [51 For low Coulomb 
interaction (see plot in solid thick line) the induced ef- 
fective pumping potential is small and the behavior of 
the dc-current is qualitatively the same as the one corre- 
sponding to the non-interacting case. For increasing U, 
the induced pumping amplitude Veff becomes sizable. 
Within the self-consistent procedure, this is a complex 
function of tp with real and imaginary parts V^rr and 

V^fj displaying the functional structure of ip suggested 
by the perturbative solution leading to Eq. (|^^ . The en- 
suing curve J'^'^ vs (p also shows the pattern predicted by 
Eq. (|22p . Notice that the dc-current, as well as Veff can 
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effects, which become relevant at all frequencies at finite 
T. The evolution of the dc-current as a function of Qq 
as the temperature grows within a range < T < Tk 
is analyzed in Fig. [71 At first glance, it becomes ap- 
parent how the dominant contribution at low fio is J'^^ 
at the lowest T, while it turns to be J*" for the highest 
ones. The elastic contribution vanishes as the Kondo res- 
onance disappears while the inelastic one increases due 
to the corresponding grow of F*" . The latter increment 
translates into the inverse time t^^, thus increasing the 
range of flo where the behavior J'^'^ oc SIq holds. In fact, 
notice that the range of low fig where the plots of Fig. 3] 
look horizontal increases with increasing T. 



FIG. 7: (Color online) J"' /Qo (top left), J'" /ilo (bottom left) 
and J'^'^/fto (right) as functions of the pumping frequency 
r^o for 20 values of equally spaced temperatures in the range 
0.04 < T/Tk < 0.8 and Coulomb interaction U = 10. The 
lowest and highest temperatures are indicated in thick solid 
and dashed lines, respectively. The local density of states of 
the dot, po{i^), is depicted in the inset for the lowest and 
highest temperatures. Other details are as in Fig. [2] 



be fitted with an excellent degree of accuracy by func- 
tions of (f with the structure suggested by this equation 
(see symbols in Fig. |6]). An striking feature observed 
in this figure as well as in the analytical expression ((22|) 
is the breaking of the symmetry ip — > —ip in the be- 
havior of the dc-current. On general physical grounds, 
the symmetry of the problem in the case of identical 
pumping amplitudes at the two barriers indicates that 
J'^'^{(p) -^ —J'^'^(—(p). However, the many-body treat- 
ment adopted in the present work breaks such symmetry 
for high enough induced Vg//, even at the level of the 
simple self-consistent TDHF approximation, i.e. even 
disregarding self-energy effects. The self-consistent eval- 
uation of higher harmonics (we recall that we truncate 
at |fc| = 2), and additional self-energy and vertex correc- 
tions that we have not taken into account are expected to 
restore such symmetry. In any case, our results should be 
interpreted as a piece of evidence on the departure from 
the J'^'^ oc sm{ip) behavior induced by the interactions. 

So far, we have focused in the case of temperature 
T = 0, where inelastic effects play an insignificant role at 
low Qq. To finalize, we analyze in what follows inelastic 



V. SUMMARY AND CONCLUSIONS 

To conclude, we have analyzed a simple model for an 
interacting quantum pump by means of non-equilibrium 
Green's function techniques and within a second order 
self-energy approximation. We have shown that the ef- 
fective time-dependent scattering center induced by the 
interactions generates interference effects, which should 
be detected in a J'^'^ vs (p experimental curve, following 
the pattern predicted by Eq. (P^ . We have shown that 
the Kondo effect manifests itself in the J'''^ vs fio behav- 
ior, which could be also detected in future experiments. 
Below the Kondo temperature, Tr-, the Kondo resonance 
enables the elastic transport of electrons, however the 
frequency range within which J'^'^ behaves linear in flo 
is extremely narrow. As the temperature grows, inelas- 
tic scattering becomes dominant and this range becomes 
wider. 
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